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Abstract 

Ph, We present analytical results for the leading contributions of the top quark to the 

Q_ij electroweak p parameter at order GpM® and G 2 F Mfa s . The Higgs boson and the 

gauge bosons are taken to be massless in this limit. The correction of order G^pMf is 
found to be sizeable in comparison to the the leading two-loop 0{G 2 F Mf) correction, 
however it is much smaller than the subleading 0{G\M 2 M'^) correction. 



1 Introduction 



The standard model of elementary particle physics has been in agreement with all exper- 
iments performed till date. The one ingredient that has not been found yet is the Higgs 
boson. Direct searches at LEP200 have been able to put a lower limit of 113 GeV on its 
mass As the presence of the Higgs boson is necessary for the renormalizability of the 
theory, radiative processes calculated without taking into account the graphs containing 
virtual Higgs bosons would be infinite. Therefore radiative corrections should grow with 
the Higgs-boson mass. It is a well-known effect, sometimes called Veltman's theorem [0, 
that such corrections to low-energy parameters grow logarithmically with the Higgs mass. 
Therefore precision tests as performed at LEP can use the existing data to put bounds on 
the Higgs mass. An analysis of the data indicates that the Higgs boson should be light 
(< 170 GeV) ||. Therefore the Higgs boson should be found at the next generation of 

■^Permanent address: Institute for Nuclear Research, Russian Academy of Sciences, 60th October An- 
niversary Prospect 7a, Moscow 117312, Russia. 



colliders. However there is some lingering doubt as to the unavoidability of this conclusion. 
On the one hand there is a slight discrepancy between the leptonic and the hadronic data 
in the precision measurements. Taking the leptonic data only, the preferred value of the 
Higgs boson mass would be so light as to almost rule out the standard model already. On 
the other hand one should not change the model too much, as one can get rather large 
effects quickly, for instance a fourth generation appears to be definitely ruled out. The one 
option that is not quite under theoretical control is the case of strong interactions, which 
corresponds to a naively very heavy Higgs boson (Mh > 1 TeV). In this case higher order 
corrections could be large and cancel the leading order ln(M#) effects, ultimately giving 
similar corrections as a light Higgs boson. Calculations at the two-loop order have been 
performed and actually give a negligeable effect ||. The reason for this has been traced 
back to some apparently accidental cancellations between the coefficients in the Higgs mass 
dependent terms ||. These cancellations appear to be typical for the low-energy vector 
boson-propagator and triple- vector boson interactions only jJJ. Corrections to e.g. the four- 
vector boson interactions |J or the Higgs-propagator do not exhibit such cancellations. 
Therefore it may be possible that the full effect of the strong interactions shows up only 
at the three-loop level for the precision parameters at LEP. It is therefore imperative to 
push the calculations towards the three-loop order. 

Of course, to determine the Higgs mass effects it is necessary to also have the effects of 
the other particles under very precise control, as these may be dominating. A particularly 
important parameter here is the so-called p parameter which measures the relative strengths 
of the neutral and charged currents. At the tree-level this parameter is one, due to a residual 
SU(2), so-called custodial symmetry. This symmetry gets violated by the hyper-charge 
and by mass-splittings within doublets. Because the mass-splittings give contributions 
proportional to the square of the masses, the dominating effect is the contribution of 
the top quark (Ap « GpM 2 ) ||. The Higgs mass dependence is only logarithmic and 
proportional to the hypercharge coupling (Ap « g\ \u(Mh/Mw)) f|. Since the Higgs 
mass dependent effect is sub-dominant, it could be masked by higher loop effects coming 
from the top quark. Therefore the higher loop top mass dependent corrections have to 
be calculated. As the top- Yukawa coupling is not very small, one should check whether 
perturbation theory in the Yukawa coupling is actually convergent. Known corrections at 
the two-loop level are the leading (G F M^) and sub-leading (g 2 eaA ,GpM t 2 ) contributions. 
Numerically it was found that the subleading corrections are larger than the leading ones 
|i0| . This naturally brings up the question whether the three-loop effect could also be 
large. There is actually an extra enhancement of the three-loop correction, because there 
is a term behaving like n 2 c (n c is the number of colors), whereas both one and two-loop 
corrections behave like n c . If indeed the three-loop correction would be comparable to the 
two- loop one, one should consider whether there is a way to resum at least the leading order 
terms, maybe by a Borel summation. Such a resummation is not always possible, as one 
would need an alternating sign in the series. So independent of its direct phenomenological 
implications it is important to know the three-loop heavy-top correction to Ap, since it 
would give us a handle on the nature of perturbation theory in quantum field theory. Not 
very many quantities are known to such precision. 
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As far as the calculation is concerned one must remark that a complete calculation 
at the three-loop level is impossible analytically. At best this could be done numerically. 
However the situation is somewhat simplified by taking the top mass much larger than Mw 
or Mz- In this case one can perform an expansion in the external momenta. The Feynman 
integrals that have to be calculated can then be reduced to vacuum graphs. These vacuum 
graphs are easier to calculate. At the two-loop level this program can be carried out till 
the end, as the vacuum- diagrams are known for arbitrary values of the masses. This way 
the leading and sub-leading corrections due to a heavy top quark have been calculated, 
including a full Higgs mass dependence fll], HU • Numerically the Higgs mass dependence 
is quite sizeable. Unfortunately at the three-loop level vacuum diagrams are only known 
when there is one mass scale in the problem [0. That is, every line in the graph is either 
massless or has a fixed mass. Mixed Higgs-top graphs can therefore only be calculated when 
the Higgs has the same mass as the top quark, or has mass zero. Also QCD corrections can 
be taken into account analytically. In this letter we consider the case that the top mass is 
much heavier than all other masses, effectively working in the M# = limit. Within this 
limit we present the results of the calculation of the leading order M t 6 and a s M t 4 correction 
to Ap. 

This letter is organized as follows: In the next Section we state our definition of the p 
parameter and introduce our notations. Section |3] is devoted to a short discussion of the 
treatment of the diagrams and the renormalization. In Section § we present the results at 
order G F M? and G\M^a s . 



2 Definition 

The p parameter is usually defined by the ratio of the neutral and charged current coupling 
constants at zero momentum transfer: 

9 J cc (0) 1-Ap' {) 

^cc(O) is given by the Fermi coupling constant Gf determined from the p, decay rate 
whereas Jjvc(O) is measured by neutrino scattering on electrons or hadrons. This definition 
of the p parameter is principally speaking not process independent, since the radiative 
corrections depend on the value of the hypercharge of the particles in the process. However 
the leading terms in the top mass are process independent. 

The origin of the process dependence lies in the different hypercharge assignments of 
the particles that take part in the scattering cross-sections. As already mentioned in the 
introduction, the origin of a non-zero value of Ap is due to the breaking of the custodial 
SU(2) symmetry. This breaking originates from either the hypercharge coupling or from the 
doublet mass-splittings. In leading order these effects are separate, but in higher order the 
effects of the Yukawa and the gauge couplings become intricately intertwined. At the two- 
loop level there is a leading 0{Mf) correction, that is independent of the gauge-couplings 
completely. Even in the absence of the hypercharge field there is a sub leading 0(Mf) 
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contribution. Finally the exact sub-leading term becomes a function of the hypercharge 
coupling, giving rise to a more complicated formula as a function of sm 2 (9 w ) UTO] . As the 
whole correction is also strongly dependent on the Higgs mass one should be careful in 
drawing conclusions from the large size of the subleading term compared to the leading 
term. This effect may be specific to the particular observable being calculated. In particular 
there appears to be an accidental cancellation between different coefficients in the leading 
term 0(G 2 F M^) to Ap for Mjj = 0, which disappears for a larger Higgs mass and which is 
absent in the correction to Z — > bb. 

In leading order in the top mass contributions to p stem from the transversal parts of 
the (unrenormalized) self energies of the exchanged vector bosons W and Z: 

i - n^(o)/M^ o 
p i - nf*(o)/Mf o ' 1 ) 

where M Wo M Zo are the bare masses of the W and Z bosons respectively. In this form the 
results would appear infinite, but these infinities are canceled by the Higgs mass, top mass 
and W boson mass renormalizations. Corrections from vertex and box diagrams always 
involve extra powers of the weak coupling constant g wea k and are therefore suppressed by 
powers of g weak / 1 g l y V ukawa , i.e. M w /M t . 

A remark is here in order. Up to the two-loop level for the leading top mass effects the 
correction to p can be expressed in a number of equivalent ways, i.e. in terms of G\M^ 
or Mt/Myy- At the three- loop level one must be more precise. At this level the difference 
between low-energy parameters and on-shell vector boson masses plays a role even in the 
leading M t effects. In the following we will express the final answer in terms of the on-shell 
top mass and the low-energy Fermi constant Gp from //-decay 



3 Treatment of diagrams and renormalization 

There are four major contributions to the three-loop p parameter which we will discuss 
below. First, the three-loop self energy diagrams of the W and Z boson, of course. Second, 
the top quark mass renormalization at order G F Mf and GpM 2 a s . Third, the renormal- 
ization of the vacuum expectation value of the Higgs field and fourth, the Higgs mass 
renormalization. For all contributions we have limited ourselves to the leading top quark 
mass effects which simplifies the computation drastically. 

As far as the self energy diagrams of the W and Z bosons are concerned we have chosen 
to compute an asymptotic expansion w.r.t. the external momentum since we are interested 
in the zero momentum transfer limit. The treatment of the Higgs mass has been done 



using two different approaches: First, an asymptotic expansion |14[ w.r.t. the Higgs mass 
and second, setting the Higgs mass equal to zero right from the start. Both strategies lead 
to the same result as no Higgs mass dependence is present in the final results. 

Throughout this computation we have used several computer programs and packages. 



The diagrams were generated using QGRAF |L5| and the asymptotic expansion in limit of 



the small Higgs mass was performed by EXP [fB|. The resulting three-loop integrals were 
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(a) 



H 
(b) 



Figure 1: (a) A diagram for which the traces with an odd number of 75 matrixes do 
contribute; (b) a diagram which can be safely treated with the completely anticommutating 

75- 



evaluated with the help of the program packages MINCER (I7j and MATAD [y| written in 
FORM |1{J . The two- loop top quark mass renormalization constant requires the calculation 
of two- loop on-shell diagrams. We have used the packages TARCER [EDJ and 0NSHELL2 [2I| 
for the computation of these types of diagrams. 

The renormalization of the vacuum expectation value of the Higgs field is fixed by the 
requirement of vanishing Higgs tadpoles. It also influences the Higgs mass renormalization 
and the mass renormalization of the Goldstone bosons. We need this renormalization 
constant up to two-loops and the computation of the tadpoles is straightforward. The 
mass renormalization of the Higgs field is required only up to one-loop order. 

We have used the prescription of a completely anticommutating 75 in the evaluation 
of all traces. In particular, all traces with an odd number of 7s's were set to zero. This 
treatment is completely justified if and only if every such diagram vanishes on kinematical 
grounds. For instance, obviously no problem could result from 2-point amplitudes with 
only one fermion loop: the number of available independent external momenta (one) and 
Lorentz indices (two if both legs are vector bosons) is not enough to form a completely 
ntisymmetric tensor e^ u \ p . The situation can be different only for singlet type 3- loop 
diagrams comprising two fermion traces (see Fig. 1). 

Indeed, the trace occurring in one-loop vertex diagram with three vector external legs 
(Fig. 1(a)) can only be computed by using the 't Hooft-Veltman prescription for 75. For- 
tunately, in our case the only dangerous diagrams are those of Fig. 1(b). They contain at 
worst triangle subdiagrams comprising a trace of an odd number of 7s's and depending on 
only three (not four!) independent external Lorentz vectors and indices (two independent 
external momenta and at maximum one index from a vector boson leg). 

We have also successfully checked the Ward-Takahashi identities which relate the self- 
energies of the gauge bosons to the non-diagonal self-energies of the gauge-to-Goldstone 
and the Goldstone boson self-energies. 
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4 Results 



We have reproduced the one- and two-loop results for the p parameter at orders GVM 2 , 
GpMt and G F M 2 a s that can be found in IRL fFl 



Ap (2,a s M?) 



n c X t , 

n c X 2 (19 - 12( 2 ) 



a 



n c X t ^C F (---C. 



71 
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X t is defined by (M t is the top quark pole mass) 



X, 



G F M? 



(3) 
(4) 

(5) 



(6) 



8\/27r 2 

and n c represents the number of colours. ( n is used for the values of the Riemann (^-function 
at n, C(n), giving e.g. ( 2 = tt 2 /6. Cp is equal to 4/3 for the SU(3) colour group. 

For the light Higgs boson the subleading two-loop 0(G 2 F M 2 M^) corrections are much 
larger than the leading two-loop 0{G 2 F Mf) term, which is actually suppressed by accidental 
cancellations [flO | . So they should necessarily be taken into account. In contrast to the 



leading corrections the subleading M 2 two-loop corrections depend on the processes used to 
define the neutral and charged currents in Eq. (|l]). If the charged current coupling constant 
is defined through the //-decay and the neutral current coupling constant is defined through 
the v^e scattering, then the subleading two-loop corrections in the massless Higgs limit are 
given by 



Ap 



(2,M?) 



nX 2 ^ 
UcAt M 2 



+ 7 + 



11 3 319 s 2 



2 
6 



6c 2 h + ir 2 
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4r \n(dr) + 21 - 16 r) In 



56 s< 



(7) 



with I 3 representing the isospin of the electron target (J 3 = —1) [fLOR . s and c are abbre- 
viations for the MS quantities sin^ and cos 6^, respectively, given by c 2 ~ ^(1- n c X t ) 
with the on-shell parameter c 2 = M^/M§ . In the limit of vanishing hypercharge (s = 0) 
the corresponding correction looks like 

A p(2M?) 



s=0 



-1 



7vr 2 



+ 6/3 + 2I In 



M? 



(8) 



Numerically the subleading corrections (0-|D are about 20 times larger than the leading 
two-loop correction (^) and are about 4% of the one-loop correction (Q). 
The three-loop results read 

Ap (3,M t 6 ) = 

Xfn c (68 + 729S 2 + 36D 3 + 96C 2 In 2 + 6C2 - 612C 3 + 324( 4 - 725 4 ) 

6572 4374 n 1472 AAn \ 

S 2 + —(2 + 440(3 ) (9) 
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xf 


X 2 Ml 


X t a s 


XI 


Xla s 


X t a 2 s 


9.6 • 10" 3 


-2.3 • 1(T 5 


-4.1 • 10~ 4 


-9.5 • 10" 4 


8.2 • 1(T 6 


1.0 • 10" 6 


-1.7- 10- 4 



Table 1: Size of different contributions to the p parameter according to Eq. (|l^). The 
numbers are obtained with a s (M t ) = 0.109. 



and 



Ap 



(3,M t 4 a s ) 



7T 



185 729 



151 



Xfn c —C F — + — S 2 - 48( 2 In 2 - —( 2 + 29C 3 - 24( 4 + 125 4 . (10) 
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The constants S 2 , -D 3 , and B4 are given by (see, e.g. JT8 

4 



S 2 



Cl 2 



9^ 
6C3 " ^C4 - 6 



0.260434 . . . 

2 



Cl 2 



7T 



16Li 4 



3.027009 . . . 
13 



4C 2 ln 2 2 + -ln 4 2-yC4 



-1.762800... 



with Cl2(7r/3) = Im[Li 2 (e l7r / 3 )]. Li 2 and Li 4 represent the Di- and Quadrilogarithm, respec- 
tively. 

Numerically, we find 



Ap 



n c X t - 0.7392n c X 2 t - X 2 t ^§ n c (30.029 - 21 In 

+ n c (10.1466 + 24.3669 n c ) Xf 

+ —n c (-2.8599 X t + 0.9798 X 2 ) +■■■ 



'Ml 
MI 



(12) 



and, for n r 



Ap 



M 2 ( 

3 X t - 2.2176 X 2 -X 2 ^[ 90.1 - 63 In 

Mf \ 



+ 



f 

a i 

7T 

43.782 



'Mz 
MI, 



249.74 X? 



8.5797X, + 2.9394 X; 



X t + 



(13) 



Here the subleading two-loop correction is given again for the process of u^e scattering 
and in order to see the pure effect of the top quark mass on Ap the value s = was 
taken. For the realistic value of s 2 = 0.2311 the subleading two-loop corrections are 
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-XfM|/M t 2 (78.7-51.9 ln(M|/M t 2 )). For completeness we have also added to the r.h.s. of 
13|) the well known term of order GpMfa 2 @ . 



The size of the various terms in Eq. (|13|) is shown in Table 1. A few remarks are in order. 
First, the coefficients of the perturbative series are alternating for both pure electroweak 
and mixed electroweak/strong corrections. Second, the size of the three-loop electroweak 
coefficient is quite large and results in a correction to the leading two-loop contribution of 
approximately 36 % (for X t w 3.2 ■ 10~ 3 ). The size of the three-loop correction relative 
to the much larger sub leading two-loop M f 2 corrections (0-§) is however about 2%. The 
mixed 0(G 2 F Mfa s ) term is numerically less than the purely electroweak 0(G F Mf) one by 
almost an order of magnitude. Finally, the major part of the 0(Xf) contribution to the p 
parameter is given by the leading term in the large n c approximation. 

With respect to the question of the convergence of perturbation theory in the Yukawa- 
coupling X t the results give only limited information. Because of the dependence on n c 
one can actually study two different limits, selecting different subsets of graphs. The first 
possibility is to take the graphs with the largest number of fermion- loops into account, 
which is one less than the number of loops, starting at the two-loop level. Formally 
one can project out these graphs by taking X t « l/n c and take the limit n c — > oo. 
With the exception of one-loop all other loops are of the same order. One can thus 
only meaningfully compare the two- loop with the three- loop level. These become equal 
for n c X t ~ 0.03, corresponding to M t ~ 300 GeV. This would naively imply an early 
breakdown of perturbation theory. Such a conclusion however would be premature, as it is 
known that the two-loop p parameter has a large cancellation between the different terms, 
that have not been observed for other quantities. 

The other limit is taking into account only graphs with one fermion-loop, which one can 
describe formally by taking the limit n c — > 0. Here one has three terms in the perturbation 
expansion available and can therefore form the [l,l]-Pade approximant. This gives 

v 1 + 12.987X* . , 

Ap = n c X t . 14 

p l + 13.7265X t v ; 

For a large top mass (Xt — > oo) this would imply a saturation of the radiative correction 

Ap = 0.946 A Pl _ loop . (15) 

However also this early saturation could be an artefact, coming from the anomalously small 
two-loop leading term. All in all one gains little insight here without having a way to sum 
high-loop graphs. Also a comparison with other quantities would be necessary as well as 
a study of the Higgs mass dependence. 

Conclusion 

We have analytically computed the leading contributions of the top quark to the elec- 
troweak p parameter at order GpMf and GpM 2 a s . The correction of order G F Mf is 
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sizeable in comparison to the leading two-loop correction. As we have mentioned in the 
Introduction the three-loop leading correction in principle could have happened to be of 
the same order of magnitude as the subleading G\M^M\ two-loop correction. We have 
found, however, that in fact it is still much smaller than the subleading two-loop one. 
Thereby we have significantly reduced the theoretical uncertainty of the results of elec- 
troweak higher order calculations coming from the enhanced three-loop large top quark 
mass contributions. 
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